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Abstract 

Vacuum expectation value of the energy-momentum tensor and the vacuum interaction 
forces are evaluated for a massive scalar field with general curvature coupling parameter 
satisfying Robin boundary conditions on two codimension one parallel branes embedded 
in (D + l)-dimensional background spacetime AdSoi+i x 5] with a warped internal space 
S. The vacuum energy-momentum tensor is presented as a sum of boundary-free, single 
brane induced, and interference parts. The latter is finite everywhere including the points 
on the branes and is exponentially small for large interbrane distances. Unlike to the purely 
AdS bulk, the part induced by a single brane, in addition to the distance from the brane, 
depends also on the position of the brane in the bulk. The asymptotic behavior of this part is 
investigated for the points near the brane and for the position of the brane close to the AdS 
horizon and AdS boundary. The contribution of Kaluza-Klein modes along E is discussed 
in various limiting cases. The vacuum forces acting on the branes are presented as a sum of 
the self-action and interaction terms. The first one contains well known surface divergences 
and needs a further renormalization. The interaction forces between the branes are finite 
for all nonzero interbrane distances and are investigated as functions of the brane positions 
and the length scale of the internal space. We show that there is a region in the space of 
parameters in which these forces are repulsive for small distances and attractive for large 
distances. As an example the case S = S^^ is considered. An application to the higher 
dimensional generalization of the Randall- Sundrum brane model with arbitrary mass terms 
on the branes is discussed. Taking the limit with infinite curvature radius for the AdS bulk, 
from the general formulae we derive the results for two parallel Robin plates on background 
of X E spacetime. 



PACS numbers: 04.62. -Fv, ll.lO.Kk, 04.50.-hh 



1 Introduction 

Recent proposals of large extra dimensions use the concept of brane as a sub-manifold embed- 
ded in a higher dimensional spacetime, on which the Standard Model particles are confined. 
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Braneworlds naturally appear in string/M-theory context and provide a novel setting for dis- 
cussing phenomenological and cosmological issues related to extra dimensions. The model in- 
troduced by Randall and Sundrum W is particularly attractive. The corresponding background 
solution consists of two parallel flat 3-branes, one with positive tension and another with neg- 
ative tension embedded in a five dimensional AdS bulk. The fifth coordinate is compactified 
on /Z2, and the branes are on the two fixed points. It is assumed that all matter fields are 
confined on the branes and only the gravity propagates freely in the five dimensional bulk. In 
this model, the hierarchy problem is solved if the distance between the branes is about 37 times 
the AdS radius and we live on the negative tension brane. More recently, alternatives to confin- 
ing particles on the brane have been investigated and scenarios with additional bulk fields have 
been considered. 

From the point of view of embedding the Randall-Sundrum model into a more fundamental 
theory, such as string/M-theory, one may expect that a more complete version of this scenario 
must admit the presence of additional extra dimensions compactified on an internal manifold. 
From a phenomenological point of view, higher dimensional theories with curved internal man- 
ifolds offer a richer geometrical and topological structure. The consideration of more general 
spacetimes may provide interesting extensions of the Randall-Sundrum mechanism for the ge- 
ometric origin of the hierarchy. Spacetimes with more than one extra dimension can allow for 
solutions with more appealing features, particularly in spacetimes where the curvature of the 
internal space is nonzero. More extra dimensions also relax the fine-tunings of the fundamental 
parameters. These models can provide a framework in the context of which the stabilization 
of the radion field naturally takes place. In addition, a richer topological structure of the field 
configuration in transverse space provides the possibility of more realistic spectrum of chiral 
fermions localized on the brane. Several variants of the Randall-Sundrum scenario involving 
cosmic strings and other global defects of various codimensions have been investigated in higher 
dimensions (see, for instance, P]-jl4j and references therein). In particular, much work has been 
devoted to warped geometries in six dimensions. 

Motivated by the problems of the radion stabilization and the generation of cosmological 
constant, the role of quantum effects in braneworlds has attracted great deal of attention (see, 
for instance, references given in Iii this paper we continue the investigation of the local 

quantum effects induced by two codimension one parallel branes embedded in the background 
spacetime AdSo^+i x S with a warped internal space S (for the investigation of the vacuum 
energy-momentum tensor in the geometry with AdSo+i bulk see Refs. |16( I17j). The quantum 
effective potential and the problem of moduli stabilization in the orbifolded version of this model 
with zero mass parameters on the branes are discussed recently in Ref. |18j . In particular, it 
has been shown that one loop-effects induced by bulk scalar fields generate a suitable effective 
potential which can stabilize the hierarchy without fine tuning. In the previous paper we have 
studied the Wightman function and the vacuum expectation value of the field square for a scalar 
field with an arbitrary curvature coupling parameter obeying Robin boundary conditions on the 
branes. For an arbitrary internal space S, the application of the generalized Abel-Plana formula 
allowed us to extract form the vacuum expectation values the part due to the bulk without 
branes and to present the brane induced parts in terms of exponentially convergent integrals for 
the points away the branes. An application of the general results to the model with S = 
was discussed. Here we consider the vacuum expectation values of the energy-momentum tensor 
and the vacuum interaction forces between the branes. In addition to describing the physical 
structure of the quantum field at a given point, the energy- momentum tensor acts as the source 
in the Einstein equations and therefore plays an important role in modelling a self-consistent 
dynamics involving the gravitational field. 

The paper is organized as follows. In the next section, by using the expression for the Wight- 
man function from ' 15J , we derive formula for the expectation values of the energy-momentum 



2 



tensor for a single brane geometry in both regions on the left and on the right of brane. Various 
limiting cases are discussed in the general case of the internal manifold. The corresponding 
quantities for the geometry of two branes are investigated in section |31 The case of different 
Robin coefficients on separate branes is considered and an application to the higher dimensional 
version of the Randall-Sundrum braneworld is discussed. In section |1] we study the vacuum 
interaction forces between the branes as functions on brane positions and AdS radius. The 
application of the general formulae to the example with the internal space T, = with the 
further generalization to the case S^^ is discussed in section |SJ Section |S1 contains a summary 
of the work and some suggestions for further research. 

2 Casimir densities for a single brane 

In this paper we are intereseted in local quantum effects arising from a scalar field (p{x) prop- 
agating on background of {D + l)-dimensional spacetime with topology AdS^i+i x S and the 
line element (we adopt the conventions of Ref. 19 for the metric signature and the curvature 
tensor) 

ds'^ = QMNdx^Ux^ = e-^'^^yrif.^dx^dx'' - e-^'^^y-fikdX'dX'' - dy'^, (2.1) 

where 77^0- = diag(l, —1, . . . , —1) is the metric tensor for Di -dimensional Minkowski spacetime 
R{Di-hi) ^ j.^ is the inverse AdS radius, and the coordinates X*, i = 1, . . . , D2, cover the manifold 
T,, D = Di + D2. Here and below the upper-case latin indices M,N range from to D and 
fj,, a = 0,1, Di — 1. In addition to the radial coordinate y we will also use the coordinate 
z = e^'^y /kj:), in terms of which line element (|2.1I) is written in the form conformally related to 
the metric in the direct product spacetime R^^^'^^ x S by the conformal factor {kDz)~'^ . The 
solutions to Einstein equations of type (|2.1() have been considered in Refs. [511^1 E]- For the 
scalar field with curvature coupling parameter C, the field equation has the form 

{V^^Vu + m^ + CR)v{x) =0, (2.2) 

with R being the scalar curvature for the background spacetime and Vm is the covariant deriva- 
tive operator associated with the metric tensor qmn- For minimally and conformally coupled 
scalars one has C = and Q = C,o = {D — 1)/ (41)), respectively. In the case of the bulk under con- 
sideration, the corresponding Ricci scalar is given by formula R = -D{D + l)kl - e^'^^y Rf^^-^, 
where -R(^) is the scalar curvature for the metric tensor 7jfc. In the discussion below we will 
assume that the field obeys Robin boundary conditions 

(4 + Bydy^ ip{x) = 0, y = a,b, (2.3) 

with constant coefficients Ay, By, on two parallel branes of codimension one, located at y = a 
and y = b, a < b. The expressions for these coefficients in a higher dimensional generalization 
of the Randall-Sundrum brane model will be given below. The z-coordinates of the branes 
we denote by zj = e^^^ /kf), j = a,b. Robin type conditions are an extension of Dirichlet and 
Neumann boundary conditions and appear in a variety of situations, including the considerations 
of vacuum effects for a confined charged scalar field in external fields, spinor and gauge field 
theories, quantum gravity and supergravity. These boundary conditions naturally arise for 
scalar and fermion bulk fields in braneworld models. The Wightman function and the vacuum 
expectation value (VEV) of the field square for the geometry of two branes in the bulk with line 
element ()2.1() are investigated in our previous paper Here we are interested in the VEV of 
the energy-momentum tensor and the vacuum interaction forces between the branes. 
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The VEV of the energy-momentum tensor can be evaluated by substituting the Wightman 
function and the VEV of the field square into the formula 



(0|rMiv(x)|0) = lim dMd'j^{0\ip{xMx')\0) 

1 



+ 



C - ^ ) Sa/jvVlV^ - CVa/Vjv - CRmn 



{^W\xm, (2.4) 



with the components of the Ricci tensor (here and below the components of tensors are given 
in the coordinate system with the radial coordinate y) 

Rfia = -Dkjjg^^, Rik = R(^)ik - Dk'jjgik, Rdd = Dkj), (2.5) 

where R(^-)ik is the Ricci tensor for the metric 7,^. This corresponds to the point-splitting 
regularization technique for the VEVs. First let us consider the geometry of a single brane 
located at y = a. In the region y > a the corresponding Wightman function is determined by 
the formula ^3] 

13 



duu — Ju{uz)Jy{uz') / duu ^, , — 



(t_t')^^2_fc2j I (2.6) 



K^(uz)K^{uz') , 
^=^= — cosh 

V? — k'^ 

with kj^ = ^ k'^ + A^, A; = |k|, and k is the wave vector in the subspace with spatial 

coordinates x = {x^ , . . . ,x^i~^). In formula (|2.6j) . Ju{x) is the Bessel function, Iu{x), Ky{x) are 
the Bessel modified functions with the order 



L>2/4-L>(D + l)C + m2/A;|), (2.7) 

'4)i3{X) are the eigenfunctions for the operator A(^) + C-^(7) with eigenvalues — A^, and A(^) is 
the Laplace-Beltrami operator for the metric 7jfc. Here and below for a given function F{x) we 
use the notation 

F^^\x) = AjF{x) + BjxF'{x), (2.8) 
with the coefficients related to the constants in the boundary conditions by the formulae 

Aj = Aj + BjkDD/2, Bj = Bjko. (2.9) 

In formula 1)2. 6|) , the part with the first integral in the figure braces is the Wightman function 
for AdS Di+i X 5] spacetime without branes, and the part with the second integral is induced 
in the region y > a by the presence of the brane. The formula for the Whightman function in 
the region y < a is obtained from Eq. (|2.6|) by the interchange of the Bessel modified functions, 

«^ Ki,. Substituting the Wightman function and the expression for the VEV of the field 
square from Ref. into formula (|2.4)) . the VEV of the energy-momentum tensor is presented 
in the form 

(0|Ti^|0) = (Ti^)(0) + (Ti^)('^), (2.10) 

where 

(ri^)W = ^i^r(i-^)5]|^,(x)|2/ duuiu^ + xl)^-^Fl,-}^[Muz)], (2.11) 

(47r)~ V ^ / Jo 
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is the VEV for the energy-momentum tensor in the background without branes, and the term 



(ril)(")= J]|#Wl'(rM)?^ (2.12) 

is induced by a single brane at z = Za, with the contribution of the given KK mode (3 along E 
determined by the formula 



/rpN\{a) 
\-'-Mlp 



'Di 



, . 2 .2^^^-l A"'\uZa) (+)N . , 

duu{u -A^)2 ^ F^j^; {Ky(uz)\. 

Kl '(uZa) 



To derive this formula we have used the relation 

uk^'-^jiu) _ r(§)r(i-g) 
2r (§ + 1 - 



(2.13) 



dk 



du- 



(n2 - )-f ' 



(2.14) 



with s = —1/2,1/2 and n = Di — 1. For a given function g{v), the functions -F^^^[ff(^^)] in 
formulae ()2.11|) and 1)2. 13|) are defined by the relations 



1 



d 



'^M 7 - C z'g'{v)vpiX) + 2v—F[g{v)] + 



dv 
v'g'\v) 



DiiC-l/A) 



g\v) ) ,(2.15) 



+D{AC - l)vg{v)g'{v) + 



2m^ 



u'±v'\g\v) 



for the components in the AdS part, and by the relations 



z?(±)'=r ^ M 2 2/ ^ , 1 



d 



+ -5f{l-A0v—F[g{v)\ 



(2.16) 

(2.17) 
(2.18) 



with t\^{X) = —^^Hi3ii{X), for the components having indices in the internal space. In these 



expressions we use the following notations 

F[g{v)] = vg{v)g'{v) + \{D + ^-^g\v), 

MX) - ' vp{x)--^ IMX)^' 



tpik{x) = V(^),^;3(x)V(^)fcV';(^) + 



ik 



\MX)\\ 



(2.19) 
(2.20) 

(2.21) 



where V(^)j is the covariant derivative operator associated with the metric tensor 7jfc. By using 
the equation for the Bessel modified functions, the derivative in expressions 1)2. ISp and H2.18p 
can also be presented in the form 



v^F[g[v)]=v^g'''{v)+[D + 



4C 



4C-1 



vg{v)g'{v) + {v +v )g {v) 



(2.22) 



for g{v) = Iy{v),K^{v). 
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By a similar way, for the VEV induced by a single brane in the region z < Za one obtains 



(ril)Jr^ = - \ rduuiu^ - A|)^-i^^i^F(tf [/.M]. (2.23) 

^ (47r)^r(^)A^ ir{uza) ^ 

In the case of the purely AdS bulk without an internal space, from formulae (|2.13p and 1)2. 23() 
we obtain the results derived in Ref. J7j- we see from the formulae given above, for the 
general case of the internal subspace S, the vacuum energy-momentum tensor is non-diagonal. 
In the case of homogeneous internal space the contributions into the energy-momentum tensor 
coming from the terms containing ?7/3(X) and rfp{X) vanish. In particular, in this case one has 

(ri))(") = 0. Note that, unhke to the AdSo+i bulk, VEVs fTT^ and (fT^ in addition to 
the ratio z/ Za depend also on the absolute position of the brane. For a one-parameter internal 
space of linear size L one has ^ 1/L. In this case the brane induced part in the vacuum 
energy- momentum tensor is a function on the ratios L/ Za and z/ Za- The latter is related to the 
proper distance of the observation point from the brane by the equation 

z/za = e^°^y-''\ (2.24) 

From the point of view of an observer residing on the brane, the physical size of the subspace S is 
La = Le~^°"' and the corresponding KK masses are rescaled by the warp factor, A^"^ = X^e^^"'. 
Now we see that the brane induced part in the vacuum energy-momentum tensor is a function 
of the proper distance from the brane and on the ratio La/il/kn) of the physical size of the 
internal space S for an observer living on the brane to the AdS curvature radius. For a fixed 
value of the distance of the observation point from the brane, the scaling of L is equivalent to 
the shift of the brane position. 

Note that if we consider a quantum scalar field with mass m on background of spacetime 
i? X S with the line element ds"^ = dt^ — ^i^dx^dx^ , when the VEVs of the spatial components 
of the energy-momentum tensor are expressed through the tensor (|2.21j) by the formula 

(Oe|T.,|0^) = 1^ Yf^^ , (2.25) 

where |0s) is the amplitude for the corresponding vacuum state. The expression for the brane- 
free part 1)2. llj) is divergent and needs regularization with further renormalization. As it was 
discussed in Ref. J3] for the case of the field square, this can be done combining the zeta function 
technique for the series over (3 and dimensional regularization for the integral over u. The series 
over [5 is expressed in terms of the local zeta function for the operator A(^) -|- C-R(7) ~ ""^^ and 
its derivatives. In section |S1 to evaluate the brane- free part of the vacuum energy- momentum 
tensor for the example with S = S^, we will use the Abel-Plana formula which allows to extract 
the part corresponding to the bulk AdSo+i and to present the remained finite part in terms of 
exponentially convergent integrals. The brane-induced parts (|2.13|) and (|2.23|) are finite for the 
points away the brane and, hence, the renormalization procedure is needed for the boundary-free 
part only. 

For the comparison with the corresponding results in the case of bulk spacetime AdSo^+i 
when the internal space is absent, it is useful in addition to the VEVs (|2.13|) and (|2.23|) to 
consider the VEV integrated over the subspace S: 

(rM)!lgrated = e'^^^^^^ |^ (T,^) ^ . (2.26) 

Note that due to the warp factor in this formula the volume of the extra space S exponentially 
decreases as one moves toward the AdS horizon. Comparing this integrated VEV with the 
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corresponding formula from Ref. jl7j , we see that for a homogeneous subspace S the contribution 
of the zero KK mode (A^ = 0) into the integrated ^- and ^-components of the brane induced 
vacuum energy- momentum tensor differs from the corresponding formulae in the AcISd-^^^i bulk 
by the replacement Z) — > Di in expressions H2.15() . ()2.16|) and in the definition of the order for 
the modified Bessel functions: for the latter case — > i^i with vi defined by Eq. 1)2 .7(1 with D 
replaced by Di. Note that for a scalar field with ^ < (^d^u^^^i one has i' > i^i- In particular, 
this is the case for minimally and conformally coupled scalars. 

It can be checked that the both boundary-free and brane induced parts in the VEV of the 
energy-momentum tensor obey the continuity equation V^T^ = 0, which for the geometry 
under consideration takes the form 



d_ 

dz 
d_ 

dz 



-DrrD 



T^) + DiT^ + TI + — V 
1 



Tii 



0, 



+ 



D 



0, 



(2.27) 
(2.28) 



with V(^)jr|) 



di{^/^Tl)) / . In particular, for a homogeneous internal space the second 
equation is satisfied trivially and the last term on the left of the first equation vanishes. Note 
that we also have the relation V(^)fet^j = 0. By using the equation for the eigenfunctions ipjj^X), 
it can be seen that the following trace relation takes place 



-{D2 - l)CD,-ir?/3(X) - A 



(2.29) 



On the base of this relation we easily verify that for a conformally coupled massless scalar the 
brane induced VEVs (|2.13jl and (|2.23j) are traceless. The trace anomalies are contained in the 
boundary-free part only. 

To clarify the dependence of the boundary induced part in the VEV of the energy-momentum 
tensor for general case of the internal space, it is useful to consider various limiting cases when 
the corresponding formulae are simplified. In the discussion below for these cases it is convenient 
to introduce the following functions 



c 



U 



+ (4C - l)u^5{ 



(4C-l)n^5«5' 



(2.30a) 
(2.30b) 
(2.30c) 

(2.30d) 



\Mx)r 

with / = 0, 1. Here and below we define rS-^"^ = 1 for the region y > j and n^^^ = —1 for the region 
y < j- First of all, as a partial check of the results derived above for the VEV of the energy- 
momentum tensor let us consider the limit 0. This corresponds to a single Robin plate on 
the bulk 72(^1'!) x S. For A;^, 0, from (HTj) we see that the order u of the cylindrical functions 
in formulae 1)2. 13() and 1)2. 23(1 is large. Introducing the new integration variable v = u/v, we can 
use the uniform asymptotic expansions of the Bessel modified functions for large values of the 
order (see, for instance, [20] )• To the leading order this gives: 



irpN\{a) 



(4vr) 



du (n^ 



f/3 



-Y.\Mx 

/3 

~2u\y—a 



(2.31) 
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where vp = J m? + Ao and in the expression for F^^{u, z) we should take kj:,z = 1 (as we are 



considering the hmit A;^) — > 0). In Eq. H2.31() we have introduced the notation 

^^H= / , uM^ ' J = «'^- (2-32) 



u 



In particular, for a homogeneous internal space, from (|2.3ip we see that |^-component of the 
brane induced part in the VEV of the energy- momentum tensor vanishes. For Dirichlet (Neu- 
mann) scalar with < (^jj^ the corresponding energy density is negative (positive). In the special 

case of bulk when the internal space E is absent, formula 1)2. 31(1 for the VEV 
incides with the result previously derived in Ref. |21j . 

For large values u, the integrands in formulae (|2.13j) and (|2.23|) behave as u^ig-'^^l^-^a.l 
and the integrals converge for the points away the branes. For the points on the branes the 
integrals are divergent, leading to the divergent VEVs of the energy- momentum tensor. These 
divergences are well-known in quantum field theory with boundaries and are investigated in 
general case of boundary geometry. To remove them more realistic model for the brane is 
needed (see, for instance, the discussion in Ref. and a model with the finite thickness brane 
in Ref. ). To find the asymptotic behavior of the vacuum energy- momentum tensor we note 
that near the brane the main contribution into the u-integral comes from large values u, and 
the Bessel modified functions can be replaced by their asymptotic expressions for large values of 
the argument [20]. Assuming kD\y — a\ <C 1 and Xfs\z — Za\ <C 1 (note that the second condition 

can also be written in the form A^"^|?/ — a| <C 1, the distance from the brane is much less than 
the physical length scale of the KK mode), to the leading order for the contribution of a given 
KK mode along S one finds the following result 

/T^vW ^ r l^^i + ^\ <Ba){knz,r+^F^^ 
where K{Bj) = 26oBj — 1 and the coefficients for separate components are defined by the formulae 



F^^ = D,6';,{C-Cd,), F^^b = D(C-Cd)(^--1J, (2.34a) 

4D={\-c)knZa(.z-Za)v}3{X), F^^^ = Z^^^f (^C " ^) • (2.34b) 

As the renormalized boundary-free energy-momentum tensor is finite on the brane, we conclude 
that near the brane the total vacuum energy-momentum tensor is dominated by the brane 
induced part and has opposite signs for Dirichlet and non-Dirichlet boundary conditions. Near 
the brane ^- and ^-components of this tensor have opposite signs in the regions y < a and 
y > a. For ( > (d-^ {C. < Cdi) the energy density is positive (negative) for Dirichlet boundary 
condition and is negative (positive) for non-Dirichlet boundary condition. Note that for a scalar 
field conformally coupled in D spatial dimensions one has C = Cd > Cdi ■ If we denote by pi 
and p the vacuum effective pressures in the subspace S and in the radial direction, respectively, 
when near the brane the equation of state has the form pi = —Die/{Di — 1) for a minimally 
coupled scalar and pi = Die/D2 for a conformally coupled scalar, where e is the vacuum energy 
density. In both cases one has \p/e\ ^ 1. Note that, unlike to the case of the purely AdS bulk 
without an internal space, here the vacuum energy-momentum tensor for a conformally coupled 
massless scalar field diverges on the brane. 

Next we consider the behavior of the brane induced VEV at large distances from the brane. 
For nonzero KK modes along S assuming z 3> A^^ we see that for the region of integration in 
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Eq. one has uz ^ 1 and the MacDonald function can can be replaced by its asymptotic 

expression for large values of the argument. The main contribution into the integral comes from 
values near the lower limit of integration and to the leading order the brane induced part is 
estimated by the formula 



N\{a) 
Ml 13 



/5 Ri'^XzaXp) 



(2.35) 



Hence, the contribution of the nonzero KK modes exponentially vanishes when the observation 
point tends to the AdS horizon, z ^ co. For the zero KK mode {\p = 0) assuming z » Za^ we see 
that the main contribution into the integral in Eq. 1)2. 1 3 j) comes from u for which uZa <C 1 and 
we can replace the Bessel modified functions with this argument by their asymptotic expressions 
for small values of the argument. The remained integral is evaluated by the standard formula 
for the integrals involving the square of the MacDonald function (see, for instance, [23]). For a 
homogeneous internal space to the leading order this yields to the following formula 



irpN\{a) 



with the notations 



k^+^z^^iDi + 2u) r (^ + 2u) T{R^ + u) 



2Di 



vr 



'Caif) 



2v 

^\ F, 



N 



F" 



^D2 



F, 



L»i +4i/ 

o_ DiF^ + Fi 
{Di + 2v) 



D 



where 



(±) 



n±2v 



(2.36) 

(2.37a) 
(2.37b) 

(2.38) 



In particular, for a conformally coupled massless scalar one has F£ = 0. From it follows 

that for the zero mode the brane induced VEV near the AdS horizon behaves as z^'^~'^'^ . In the 
purely AdS bulk {D2 = 0) this VEV vanishes on the horizon for u > {). For an internal spaces 
with D2 > 2i/ the VEV diverges on the horizon. Note that for a conformally coupled massless 
scalar and D2 = 1 the boundary induced VEV takes nonzero finite value on the horizon. The 
VEV from the zero mode integrated over the internal space (see Eq. ()2.26|) ) vanishes on the AdS 
horizon for all values D2 due to the additional warp factor coming from the volume element. In 
particular, the contribution of the brane into the total vacuum energy per unit surface on the 
brane in the region [z, 00), z > Za, is finite. 

For large distances from the brane in the region y < a one has z <C Za,l/Xp. The main 
contribution into the integral in Eq. (|2.23jl comes from the region u < 1/za in which uz <^ 1 
and we can replace the function I^{uz) by the corresponding expression for small values of the 
argument. As a result to the leading order we obtain the formula 



N\{a) 



M 



with the coefficients defined by the relations 

'7k _ n,,rfc jpD 



^K^%iZa)_ 
&\uZa) ' 

2^1 



F; = 21.6;, Ff = 2v5l F^ = -D, Fh = kDZ%{X). 



(2.39) 



(2.40) 
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In this limit the brane induced vacuum stresses in the directions parallel to the brane are 
isotropic. In the case of Dirichlet boundary condition the corresponding energy density is neg- 
ative for both minimally and conformally coupled scalars. As we see for fixed values ki), Za, 
Xp the brane induced VEV vanishes as 2;^+^^^ for diagonal components and as 2^+^'^"'"^ for the 
^-component when the observation point tends to the AdS boundary. In particular, the contri- 
bution of the brane to the total vacuum energy per unit surface on the brane in the region [0, z], 
z < Za, is finite (near z = the integrand in the corresponding integral over z behaves as z'^'^~^). 
The integral in Eq. (|2.39|) is simply estimated in two subcases. For the case z <C A^^ <C Za 
the main contribution comes from the lower limit and to the leading order the integral is pro- 
portional to {Xp/za)^^^^X'p^e~^'^i^^'^. As a result the brane induced vacuum energy-momentum 
tensor is suppressed by the factor (Xpz)^'^^'^^'^'^ (z/ Za)^^^'^e~'^'*'i^^'' . In the opposite limit for Za, 

2 <C Za <C A^^, which corresponds to small KK masses, A^'^'' <C /cd, the lower limit in the integral 
can be replaced by and to the leading order the contribution of the mode with a given /3 does 
not depend on A/3 and behaves as z^'^'^^ / z^^'^'^'^ . 

To see the convergence properties of the series over (3 in Eq. 1)2.12(1 , consider the contribution 
to the brane induced VEV of the energy-momentum tensor from the modes with large KK 
masses, Xf^z^Xf^Za 3> 1. Replacing the Bessel modified functions by asymptotic expansions for 
large values of the argument, to the leading order one finds 

D + l fOD „ -2u\z~Za\ 



(47r) — r (^) •''A;3 Ca[UZa) 



with the notations defined by Eq. ()2.:-{()|) . 



cJu) = ,., ^ , J = a,b, (2.42) 



and the definition for n^-'^ is given after formulae H2.3U() . Under the condition A^z^ ^ |yla/-Ba| 
or Ba = 0, we have Ca{uza) « n{Ba) and the integral can be expressed through the MacDonald 
function by using the formula 



Dl -1 



|Jd.n"(n2-A|)^-ie-^"'' = ^r(^)(^) Kn^^j2X,,), (2.43) 

where n = 0, 1. Additionally assuming A^|z — Za\ ^ 1, from ()2.41l) we find that the contribution 
of KK modes with large masses to the leading order is given by 

M (a) {kDz)^+^xf~\-^>^P\-'~''^\ 

{T§)^;^ - „ ' A,), (2-44) 

2^1+% 2 Ca{XpZa)\z - Za\ 2 

and is exponentially suppressed. In particular, the corresponding conditions are satisfied for all 
nonzero KK modes if the length scale of the internal space E is sufficiently small. In this case 
the main contribution into the vacuum energy-momentum tensor comes from the zero mode. In 
the opposite limit, when the length scale of the internal space is large, to the leading order the 
vacuum energy-momentum tensor reduces to the corresponding result for a brane in the bulk 
AdSo+i given in Ref. [T7] . 

In the limit of strong gravitational fields, corresponding to large values of the AdS energy 
scale feu when the values of the other parameters are fixed, for nonzero KK modes along S one has 
XjjZa, Xpz ^ 1. Hence, the behavior of the brane induced VEV of the energy-momentum tensor 
in this case can be estimated by formula (|2.41|1 . If in addition the condition Xji\z — Za\ 3> 1 is 
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satisfied, we have formula (|2.44l) with the exponential suppression of the brane induced part. For 
the zero KK mode the corresponding behavior of the brane induced vacuum energy-momentum 
tensor is estimated by using the asymptotic formulae for the Bessel modified functions for small 
and large values of the argument in the regions y > a and y < a respectively. To the leading 
order the components of this tensor behave like k^^'^^e^'^^'^y exp[{Di + 2u)k£){y — a)] in the 
region y < a and like k^^~^^e^^^°y eKp[2ukD{a — y)] in the region y > a. The corresponding 
quantities integrated over the internal space contain additional factor e"^^^^^ coming from the 
volume element and are exponentially small in both regions. 

Now we turn to the consideration of the limiting cases for the parameter Za determining 
the position of the brane. For small values of this parameter, Za <C 2:, for the dominant 

contribution into the li-integral in Eq. (|2.13|) one has uza <C 1 and, using the asymptotic 
formulae for the Bessel modified functions for small values of the argument, one finds 

22-Dt^2ui.D+l D2u foo „ 

Note that in this limit the distance of the observation point from the brane is large compared 
with the AdS curvature radius, koiy — CL) ^ 1, and the KK masses A^"^ are small with respect to 
the AdS energy scale, A^"'' <C /c_d. In particular, from Eq. (|2.45|) it follows that for fixed values 
ko, y, A/3 the brane induced VEV in the region z > Za vanishes as z'^'^ when the brane position 
tends to the AdS boundary, Za — > 0. Formula (|2.45l) is further simplified for two subcases. In 
the case Za <ti X^^ <^ z, or equivalently feue"*^^^^"'^^ <C A^j"'' <ti ko, the main contribution into 
the integral comes from the lower limit ant to the leading order we find 



2Di+2^7r^-Vr2(i/)ca(zy)A2 



In the limit <C ^ ^ A^ ^ (small KK masses, A^j"^ < kDC-^o^y-"^) to the leading order we can 
put in the lower limit of the integral and we obtain formula (|2.36|) . 

When the brane position tends to the AdS horizon, Za ^ 00, for massive KK modes along E 
the main contribution into the VEV of the energy-momentum tensor in the region z < Za comes 
from the lower limit of the n-integral. To the leading order we find 

{-'-M/s d] — , n /9 7\ ^ ' 



and the VEV is exponentially small. For the zero mode in the same limit, introducing a new 

(3M 



integration variable v = uZa, we expand the functions Fg^^[I,y{vz/za)] for small values of the 



argument. To the leading order this yields 



^ 2^i+2.^^r(^)r2(^ + i) \^a) Jo t\u)' 

and the suppression is power-law with respect to Za- 



/ ^-,..,Di+2u-l:^]^W 



3 Energy-momentum tensor in the region between two branes 

In the geometry of two branes we have three distinct regions of the radial coordinate: y < a, 
a < y < b, and y > b. The VEVs in the first and last regions are the same as those for a 
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single brane. In this section we will discuss the vacuum energy-momentum tensor in the region 
between two branes. Note that in the orbifolded version of the model this region is employed 
only. By using formula (|2.4|) . the VEV of the energy- momentum tensor is expressed through 
the corresponding Wightman function and the VEV of the field square investigated in Ref. jl5j . 
The Wightman function is presented in the form 

/3 



dk e^''(^-^') / duuG'^J'^ {uzj , uz)G\i^ {uzj , uz' 

QjuiuZa,UZb 







^1 



■ cosh 



k}{t-t') , (3.1) 



where kp is defined after formula (|2.6|) and j = a,b provide two equivalent representations. Here 
and in the formulae below we use the notations 

'[u)Il '{v) - Ai) '(v)Il '{u) 

and 

Gl^\u,v) = I,{v)kI^\u) - nj\u)K,{v). (3.3) 

Substituting the Wightman function into formula (|2.4j) and using relation 1)2. 14() . for the com- 
ponents of the vacuum energy-momentum tensor in the region between the branes we obtain 
the formula 

(4vr)^r(i^) V 
duu{u^ - X})^-'n,,{uZa,uzh)Fl,+f[Gi^\uz,,uz)], (3.4) 

with the functions 

^/3M^b(^)] defined by relations im^ll - lfTTHl) . where g{v) = gI^^ {uzj , v) , and 
with j = a, b. The last term on the right of this formula is finite on the brane at z = zj and 
diverges for the points on the brane z = Zj/, j' = a, b, j' ^ j. These divergences are the same as 
those for a single brane at z = Zji . As a result, if we write the VEV of the energy-momentum 
tensor in the form 

(oiri^io) = (ri^)W + (Tl^)^'^ + {tSi)^'''\ (3-5) 

j=a,h 

then the interference part (Tj^)^"'^^ is finite on both branes. In the case of one parameter 
manifold S with the size L and for a given ko, the VEV in Eq. (|3.5|) is a function on z^/za, 
L/za, and z/za- By using formula (|3.4|) and the expression for a single brane induced VEV of 
the energy-momentum tensor, we can present the contribution of the given KK mode along S 
into the interference part of Eq. (|3.5|) in the form 

(47r)^r(^) Kl >{uza)iy {uz,) _ ^ I 

K>, '(uzn) I, '(uzh) ) 
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where (T^ 



AT V (aft) . 



M//3 



is defined by the relation similar to Eq. H2.12() . and the functions Fit}^[g 



are defined by formulae (|2.15|) - H2.18() with the replacements g 
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I3M 



11,92] 



99' (5152)72, and 



g"^ — > gig2- For large values u the integrand in formula (jSISI) behaves as g-'^'^i^b-^a) 
integral is convergent for all values z including the points on the branes. For interbrane distances 
much smaller than the AdS curvature radius, {b — a) ^ l/kD, the main contribution into the 
integral in Eq. ()3.6|) comes from large values u and we can replace the Bessel modified functions 
by asymptotic expressions for large values of the argument. Additionally assuming that the 
distance is small compared with the KK length scale for an observer living on the brane y = a, 
(6- a) < 1/X^^\ it can be seen that the diagonal components of the interference part with 
M ^ D behave as [b — a)^^^^^, whereas and ^-components behave like (6 — a)~^^ . 

Now let us consider various limiting cases when the general formula for the interference part 
is simplified. First of all, in the limit fc^i — > the order of Bessel modified functions is large 
and we use the uniform asymptotic expansions for these functions. To the leading order the 
following result is obtained 



(47r 



i { — ) p Jvf3 (-a 



2\^-l 



,{ua)c},{ub)e^<^-°-) - 1 



j=a,b 



2F, 



iO)N 



0M 



{u,vp) 



kD=0 



(3.7) 



+ A^. In the expression on the right the 



with notations defined by Eq. (|2.3U|) and = 

functions F^^^ {u,v) are defined by the same relations as the functions F^'^ {u,v) from 1)2. 30() 

with the replacement n^"^ n^^'^ in formula H2.3Uc|l . and with n^^^ defined in the paragraph after 
formula 1)2. 3U() . This limit corresponds to the geometry of two parallel Robin plates in the bulk 
X S and the expression on the right of formula (|3.7|) coincides with the interference part 
(Tj(^)^(^2i,i) xs ™ VEV of the energy-momentum tensor in the region between two plates. 
This expression generalizes the corresponding result from Ref. for Robin plates in the 

Minkowski background. 

For large KK masses, ZaXp ^ 1, the arguments of the Bessel modified functions in Eq. 1)3. 6() 
are large. By using the corresponding asymptotic expansions, for the contribution of the mode 
with a given /? one finds 



irpN\{a.b) 



,D+1 



du 



Ca{uZa)cb{uZb)e'^'^('"' - 1 

{0)N, 



(4vr)^r(^ 



(3.8) 



In particular, for a one parameter internal manifold with the length scale L <^ Za, this formula 



is valid for all nonzero KK modes along S. Under the additional assumption A 



fj\z — Zj\ » 1 , the 

main contribution into the integral in this formula comes from the lower limit and we obtain 
the formula 

IBM 



lrpN\{ab) 



(47r)^ 



-2A.(^.-.)#(A„A,) 



Za) 2 



(3.9) 



Ca{>^l3Za)Cb{Xl3Zb){Zb 

In particular, it follows from here that in the expressions for the brane induced VEVs the series 
over KK modes along S are exponentially convergent. Comparing with (|2.44l) we see that for 
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large KK masses the interference part is exponentially suppressed with respect to the single 
brane parts. 

Now let us assume the conditions A^Zf, ^ 1 and ZaX/3 ^ 1. This limit corresponds to large 
interbrane distances compared with the AdS curvature radius kj^^ and is realized in braneworld 
scenarios for the solution of the hierarchy problem. In this limit the main contribution into the 
n-integral comes from the region near the lower limit and for the interference part to the leading 
order we have 



/rpN\(ab) 



-^z}\,{\pZ,)^i''\za\p) 



2D, 



vr 2 



llf\zaXp) (+)Ar 



Ki^HzaXp) 



F^l^^MK.(A^z)] , (3.10) 



for the nonzero KK modes along S. In the limit z <C A^^ <C Zf, after the replacement of the 
Bessel modified functions by their asymptotic expressions for small values of the argument, the 
corresponding formula takes the form 



/rpN\{ab) 



^-1 



m 

D 



^2 °M<^D 



+ 



Ca{l^)Cb{zbXf3)T'^{l' + 1) 
Z 



(3.11) 



where Pj^ is defined by relations 1)2. 40(1 . For A/jz » 1 and XpZa < 1, using the asymptotic 
formulae for the Bessel modified functions for large values of the argument, we obtain the 
following formula 



/rTiN\{a.b) 



{kDz)^+^X- 



2D1 



TT 2 



-\b{XpZb)z^''^ Ri'^XzaXp) 



iTi^P^ A^) - ^2c,(A^z,) ^^^^'^ 



-2X0{zi,-z) 



{1)N, 



(3.12) 



In the limit <C Zf, ^ A^^, to the leading order we can put instead of A/3 in the lower limit 
of the integral over u, and by using the asymptotic formulae for the Bessel modified functions 
for small values of the argument, it can be seen that (T^)^"^^ ~ {za/ z^Y^ g{z/ z^). In particular, 
it follows from here that the interference part in the vacuum energy-momentum tensor vanishes 
as z^^ when the left brane tends to the AdS boundary, Za —>■ 0. Under the condition z <^ Zb 
an additional suppression factor appears in the form (z/zb)^'^ for |^-component and in the form 
(z/z;,)'^!"^^"! for the other components, where ai = min(l,z^). As we see for large interbrane 
distances the interference part of the brane induced VEV of the energy-momentum tensor is 
mainly located near the brane z = Zb- 

In the higher dimensional generalization of the Randall-Sundrum braneworld based on the 
bulk AdSoi+i X S, y coordinate is compactified on an orbifold /Z2 and the orbifold fixed 
points are the locations of two L>-dimensional branes. The corresponding VEVs of the energy- 
momentum tensor are determined by the formulae given in this section with an additional factor 
1/2 and with Robin coefficients (see also [TB] ) 



Aa 
Br, 



1 



ci + ADCkj,) 



Bb 



-i-C2+AD(:kD) 



(3.13) 
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for untwisted scalar field. Here ci and C2 are surface mass parameters on the branes y = a 
and y = b, respectively. For twisted scalar Dirichlet boundary conditions are obtained. The 
one-loop effective potential and the problem of moduli stabilization in this model with zero 
mass parameters cj are discussed in Ref. jl8j . A scenario is proposed where supersymmetry 
is broken near the fundamental Planck scale, and the hierarchy between the electroweak and 
effective Planck scales is generated by a combination of redshift and large volume effects. The 
corresponding parameters satisfy the relations L < l/kn, Za ~ L, and Zb/za ^ 1. In this case 
the contribution into the interference part of the energy-momentum tensor from the nonzero 
KK modes is estimated by formula (|3.in|) . For the zero mode the interference part behaves like 
^2ukDia~b) g^j^j jg exponentially small. 



4 Interaction forces 

Now we turn to the vacuum forces acting on the branes. The corresponding effective pressure p^^^ 
acting on the brane at z = Zj is determined by ^-component of the vacuum energy-momentum 
tensor evaluated at the point of the brane location: p^^'^ = —{T£)z=Zj- For the region between 
two branes it can be presented as a sum of two terms: 



0-)=p(^Vpg,), i = a,6. (4.1) 



The first term on the right is the pressure for a single brane at z = zj when the second brane is 
absent. This term is divergent due to the surface divergences in the VEVs and needs additional 
renormalization. This can be done, for example, by applying the generalized zeta function 
technique to the corresponding mode sum. This procedure is similar to that used in Ref. JH] 
for the evaluation of the effective potential. The corresponding calculation lies on the same line 
with the evaluation of the surface Casimir densities and will be presented in the forthcoming 
paper. Here we note that the single brane term in the vacuum effective pressure will contain 
a part which depends on the renormalization scale. This part will change under the change 
of the renormalization scale and can be fixed by imposing suitable renormalization conditions 
which relates it to observables. Below we will be concentrated on the second term in the right 
of Eq. 1)4. ip . This term is the additional vacuum pressure induced by the presence of the 
second brane, and can be termed as an interaction force. It is determined by the last term on 
the right of formulae ()3.4|) evaluated at the brane location z = Zj. It is finite for all nonzero 
interbrane distances and is not changed by the regularization and renormalization procedure. In 
particular, it does not depend on the type of the regularization procedure used. All regularization 
umbiguities are involved in p^ . Substituting z = zj into the second term on the right of formula 
()3.4|) and using the relations 



dGi^\u,v) _Ai 
u 

for the interaction part of the vacuum effective pressure one has 



G\]^iu,u) = -B,, ^^^^^ (4.2) 

av v=v 



= 1 , E / duuiu" - Xl)^~%,{uZa,uzb)F^^\uz^), (4.3) 

(47r)-r(^) p Jx, 



where we have introduced the notation 

,2 



= {u' + 2^) B] - D{K - l)A,B, -A]-2[q- zp]r^p{X). (4.4) 

Below we will show that for small interbrane distances the interaction part dominates the single 
brane parts in Eq. (|4.1|) . For a Dirichlet scalar Bj = and one hasF^^'^(u) = -A^.. By using the 
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properties of the modified Bessel function it can be seen that in this case 



> and, 



hence, the vacuum interaction forces are attractive. For a given value of the AdS energy scale 
kjj and one parameter manifold S with size L, the vacuum interaction forces (|4.3j) are functions 
on the ratios z^/za and L/za- The first ratio is related to the proper distance between the branes 
and the second one is the ratio of the size of the internal space measured by an observer residing 
on the brane at y = a to the AdS curvature radius A; 7,^. Note that the term 'interaction' in the 



(7) 

discussion here and below should be understood conditionally. The quantity determines 
the force by which the scalar vacuum acts on the brane due to the modification of the spectrum 
for the zero-point fluctuations by the presence of the second brane. As the vacuum properties 
depend on the coordinate y, there is no a priori reason for the interaction terms (and also for the 
total pressures p^^^) to be equal for the branes j = a and j = b, and the corresponding forces in 
general are different even in the case of the same Robin coefficients in the boundary conditions. 

Now we turn to the limiting cases when the expressions for the interaction forces between 
the branes are simplified. First of all we consider the limit — > 0. By the way similar to that 
used before for the vacuum energy- momentum tensor, to the leading order we find 



J3) 
^(int) 



,(4vr) 



■Di 



du 



Ca{u)cb{u)e^<^-°-) - 1 



1 + 



(2C - 1/2)^1 
^2 ^oJ ^P^^) 



(4.5) 



where vp is defined after formula 1)2. The expression on the right of this formula presents the 
corresponding force acting per unit surface on the brane in the bulk geometry x S. Note 

that in this case the vacuum effective pressures are the same for both branes if the coefficients 
in the boundary conditions are the same. Moreover, for a homogeneous internal space, the 
contribution of the second term in the square brackets on the right of Eq. ()4.5() vanishes, and 
the interaction forces are the same even in the case of different Robin coefficients for separate 
branes. 

For large values of KK masses along S, XpZj ^> 1, we can replace Bessel modified functions 
by their asymptotic expansions for large values of the argument. For the contribution of a given 
KK mode to the leading order this gives 



0) 

P{mt)f} 



2{kDZj)D+^ 

■(4vr)^r(f] 



00 



du 



n2(u2-A2)^-i 



1 + 



(2C - l/2)z2 R2 



A2 - (uZjBj)'^ 



,2 



Ca{uZa)Cb{uZb)e'^^(''b 2a) 

MX) 



(4.6) 



where P^^2t)f3 determined by the relation similar to Eq. (|2.12l) . If in addition one has the 
condition Xp^Zb — Zg) ^ 1, the main contribution into the u-integral comes from the lower limit 
and we have the formula 



(i) 

^{int)/3 



)D+1 _)^i'i/2+lg-2A^(2f,-2„) 



(47r) 2 Ca{\fiZa)Cb{\pZb){zb - Za)' 



1 + 



(2C - l/2)z|i3| 
A] - (\pz,B,Y 



(4.7) 



In particular, for sufficiently small length scale of the internal space this formula is valid for all 
nonzero KK masses and the main contribution to the interaction forces comes from the zero 
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KK mode. In the opposite limit of large internal space, to the leading order we obtain the 
corresponding result for parallel branes in AdSo+i bulk |17j . 

For small interbrane distances, k£){h — a) <C 1, which is equivalent to zi^/ Za — 1 <C 1, the main 
contribution into the integral in Eq. (|4.3j) comes from large values u and to the leading order we 
obtain formula (|4.6j) . If in addition one has \p{zi) — Za) ^ 1 or equivalently X^^\b — a) <C 1, and 
assuming {h — a) <^ \Bj/Aj\ or Bj = 0, we can put in this formula = 0, and to the leading 
order one finds 

where = 1 — 2^^^ for K{Ba)K{Bh) = —1 and a^i = —1 for K{Ba)K{Biy) = 1. It follows 
from here that for small interbrane distances the interaction forces are repulsive for Dirichlet 
boundary condition on one brane and non-Dirichlet boundary condition on the another and are 
attractive for other cases. As in the limit a — > 6 the renormalized values of the single brane parts 
Pi , j = a,b, are finite, for small interbrane distances the main contribution into the vacuum 
effective pressure p^^^ comes from the interaction part. 

Now we consider the limit XpZb ^ 1 assuming that X^Za < 1. Using the asymptotic formulae 
for the Bessel modified functions containing in the argument z^, we find the following result 

2Di+ivr 2 ^ Cb{X(3Zb)Kl, [XpZa) 

This limit corresponds to the interbrane distances much larger compared with the AdS curvature 
radius and inverse KK masses, measured by an observer on the left brane: b — a ^ l/ko, 1/X^^\ 
For a single parameter manifold S with length scale L and {b — a) ^ La these conditions are 
satisfied for all nonzero KK modes. 

In the limit ZaXp ^ 1 for fixed z^Xp, by using the asymptotic formulae for the Bessel modified 
functions for small values of the argument, one finds 



^ 2-°---.g^'.f _ A^)^-#^. (4.12) 

TT 2 r [if-) vT^{u)Ca{v) lu [UZb) 



In this case the KK masses measured by an observer on the brane at y = a are much less 
than the AdS energy scale, A^"^ <C ko, and the interbrane distance is much larger than the 
AdS curvature radius. In particular, substituting A^ = 0, from these formulae we obtain the 
asymptotic behavior for the contribution of the zero mode to the interaction forces between the 
branes in the limit Za/zb ^ 1. From formulae (|4.11jl and (|4.12j) it follows that, in dependence 
of values of the coefficients in Robin boundary conditions, for large distances the interaction 
forces can be either attractive or repulsive. In particular, for Dirichlet boundary condition on 

the brane at y = b, one has ^ ^ ^ft^ i^) < ^.n^ ^(int)/3 ^ ^ l^a/^a| > and 

under these conditions for Ba 7^ we have an interesting situation when the interaction forces 
are repulsive for small distances (see Eq. ()4.8|) ') and are attractive for large distances. 

From the formulae given above it follows that in the limit when the right brane tends to the 
AdS horizon, Zb — > 00, the force P(fnt)/3 vanishes as e"^'^''^''/^^^^^ for the nonzero KK mode along 
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S and as "^^ for the zero mode. The effective pressure vanishes as z^^^^^^^^^e ^-^/'^f' 

for the nonzero KK mode and behaves as z^^"^^ for the zero mode. In the hmit when the left 
brane tends to the AdS boundary the contribution of a given KK mode vanishes as z^^'^'^ for the 
pressure ^Ifnt)/? ^ ^'a^ ^'^^ pressure V^-a^^p- -^o^ small values of the AdS curvature radius 
corresponding to strong gravitational fields, assuming XpZa ^ 1 and Xpiz}, — Zg) ^ 1, we can 
estimate the contribution of the nonzero KK modes to the vacuum interaction forces by formula 
(|4.7|) . In particular, for the case of a single parameter internal space with the length scale L, 
under the assumed conditions the length scale of the internal space measured by an observer 
on the brane at y = a is much smaller compared to the AdS curvature radius. La ^ ■ If 
La ^ one has \pZa ^ 1 and to estimate the contribution of the interference part we can use 
formulae (|4.9j) and (|4.1fl|) . and the suppression is stronger compared with the previous case. For 
the zero KK mode, under the condition A;£)(6 — a) » 1 we have Zajz^, <C 1 and to the leading 
order the corresponding interaction forces are described by relations (|4.11|) and (|4.12|) . From 
these formulae it follows that the interaction forces integrated over the internal space behave as 
fc^^"*"^ exp[(Z)i5^ + 2z/)A;£)(a — 6)] for the brane at y = j and are exponentially suppressed. Note 
that in the model without the internal space we have similar behavior with u replaced by ui and 
for a scalar field with ( < Cd+Di+i the suppression is relatively weaker. 



5 An example: S = 5^ 

To make our discussion concrete, here we consider a simple example with E = 5*^. In this 
case the bulk corresponds to the AdS^+i spacetime with one compactified dimension X. The 
corresponding normalized eigenfunctions and eigenvalues are as follows 

^^(X) = -^e^-W, Xf3 = ^\Pl /3 = 0,±1,±2,..., (5.1) 

where L is the length of the compactified dimension. The boundary-free part of the vacuum 
energy-momentum tensor can be evaluated on the base of formula 1)2. 4|) by using the Wightman 
function and the VEV of the field square from ^3]. The Wightman function is determined 
by the term in formula ()2.6() with the first integral in the figure braces. The application of 
the Abel-Plana formula to the series over /? allows to present the boundary-free part of the 
energy-momentum tensor in the form (no summation over M) 

{Ti^f^ = {T^i)%,,,+S^ / du ^S[%f!,^\u,z/L), (5.2) 
where we use the following notations 



(0) _ /I ,\_Ju-D-l)e- + u + D + l Cue- 



2„,2 



and 



(m) J (e"-l)2 e"-l 2L2 ' ^ ' 

(0) _ u ^ (n-D-l)e" + n + ZJ + l , DCue^ nV ^ ^^"^ /r, 

^(D) - 4^ (e- - 1)2 + ^> 

7r^~ r(z^ + IjF (z^ + / + ^^^j M(e" - 1) ^^L/ 
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with the hypergeometric function 1F2. The first term on the right of formula (|5.2|) is the 
corresponding tensor in AdSn+i bulk without boundaries and the second term is induced by 
the compactness of X direction. The latter is finite and the renormalization procedure is needed 
for the term {Tj^i)^^^^g^_^^ only. The renormalized value of this tensor does not depend on 
the spacetime point and is well-investigated in literature. For this reason below we will be 
concentrated on the second term. For 2: ^ L to the leading order one has 

- (Ti^)i°k,, - 26^f^-^P^CH{D + 2u)T + , (5.7) 

for M = 0, 1, . . . , -D — 1 and C-R(-2) being the Riemann zeta function. The corresponding formula 
for the component (T^^)^^^ differs from ()5.7|) by an additional coefficient —D/lu in the second 
term on the right of this formula. This directly follows from the continuity equation (j2I2ZI)- 
Hence, the second term on the right of formula 1)5. 2() vanishes at the AdS boundary as z^'^'^'^. 
For both minimally and conformally coupled cases the energy density corresponding to this term 
is positive. Note that in this limit the vacuum stresses along compactified and uncompactified 
directions are isotropic. Of course, we could expect the vanishing of the second term on the right 
hand side of Eq. (|5.2() when z/L — > 0, as this corresponds to the decompactification limit for 
the internal space. For z ^ L,hy using the asymptotic formula for the hypergeometric function 
for large values of the argument, to the leading order one finds 
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for M 7^ D—1. The corresponding formula for ^~|-component differs from ()5.8() by an additional 
coefficient —D in the second term on the right. In particular, from here it follows that the 
energy density corresponding to the second term on the right of Eq. (|5.2() is negative near 
the AdS horizon and diverges as z^ . The same is the case for (ro°)(°) as the energy density 
corresponding to the first term does not depend on z and the total energy density is dominated 
by the second term. The limit z ^ L is realized, in particular, when /c^) — > for fixed y and L. 
In this case the first term on the right of Eq. 1)5. 8() vanishes and from this formula we obtain 
the standard result for the Casimir energy-momentum tensor in 

R{D-hi) X S\ Combining the 

asymptotic formulae (|5.7j) and (|5.8|) . we see that for a scalar field with C < C^^'' (in particular, 
this is the case for minimally and conformally coupled scalars) the energy density corresponding 
to the second term on the right of formula (|5.2|) tends to zero for small values of the ratio z/L 
being positive and tends to —00 for large values of this ratio. Hence, it has a maximum at some 
intermediate value of z/L. As the first term on the right of Eq. (|5.2j) is constant, the same is 
true for the boundary-free total energy density. 

For the internal space the brane induced VEV of the energy-momentum tensor and 
the vacuum interaction forces between the branes are obtained from general formulae given in 
previous sections by the replacements 



13 f3=0 



where the prime means that the summand P = should be taken with the weight 1/2. From the 
general analysis given above it follows that under the conditions z ^ L, Za, the main contribution 
into the part of the vacuum energy-momentum tensor induced by a single brane comes from the 
zero mode and this part behaves as z^~'^'^. Comparing with 1)5. 8() . we see that for the geometry 
of a single brane, near the AdS horizon the total vacuum energy-momentum tensor is dominated 
by the second term on the right of Eq. (|5.2|1 . Near the AdS boundary, z <^ L, Za, the single brane 
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induced term in the vacuum energy- momentum tensor behaves Hke 2^+2^^ (gee Eq. (|2.39p '). By 
taking into account the estimate ()5.7() . we see that in this case the vacuum energy- momentum 
tensor is dominated by the first term on the right of (|5.2j) . For the points near the brane the 
main contribution into the energy-momentum tensor comes from the brane induced term. As 
an ihustartion, in figure Q we have plotted single brane induced VEVs for the vacuum energy 
density and ^-stress as functions on zjza for the internal space with Ljza = 1 (left panel) and 
Ljza = 2 (right panel) in the case of -D = 5 minimally coupled massless scalar field with the 
ratio of Robin coefficients Ba/Aa = 0-15. In this case for large values z/za the brane induced 
part of the energy-momentum tensor behaves as [z/za)^^- For small values of this ratio one has 
the behavior (z/za)^^- 




M =0 



M=D 
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Figure 1: Brane induced vacuum energy density (Tq)^""^ /k^~^^ and stress (Tj^)^"'^ /k^'^^ in units 
of k^^^ as functions on z/za for a minimally coupled massless scalar in D = 5. The left panel 
is for L/za = 1 and the right one is for L/za = 2. 

In figure[21we present the vacuum interaction forces in the geometry of two branes as functions 
on the size of the internal space and interbrane distance for a D = 5 minimally coupled massless 
scalar field with the Robin coefficients Ba = and Bi,/Ai, = 0.15. In this example the effective 
pressures are positive for small interbrane distances and are negative for large distances 

leading to the repulsive and attractive interaction forces respectively. 

From the point of view of embedding the Randall-Sundrum model into the string theory and 
also in the discussions of the holographic principle, the case of the internal space E = S^'^ is 
of special interest. The corresponding eigenfunctions V'/3(^) are expressed in terms of spherical 
harmonics of degree I, I = 0,1,2,.... For the internal space with radius Rq the VEVs of 
the energy-momentum tensor and the vacuum interaction forces are obtained from the general 
formulae in previous sections by the replacements 

A/3 ^ -^^/lil + D2-l)+CD2{D2-l), (5.11) 
Kq 

with the factor under the summation sign on the right in Eq. (|5.10|) being the degeneracy of 
the angular mode with a given /. 
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Figure 2: The Z? = 5 vacuum effective pressures P-^t/k^^^, j = o,,b, in units of k^~^^ for a 
massless minimally coupled scalar field as functions on L/zb and Za/zb- The values for the 
Robin coefficients are Ba = and Bb/A^ = 0.15. Left panel corresponds to j = a and right 
panel corresponds to j = b. 

6 Conclusion 

In this paper we continue the investigation of local quantum effects in higher dimensional brane 
models on the bulk of topology AcISdi+i x S with a warped internal subspace S and the line 
element (|2.1|) . The case of bulk scalar field with general curvature coupling parameter and 
satisfying Robin boundary conditions on two codimension one branes is considered. In section 
121 the model with a single brane is discussed. By using the Wightman function and the VEV 
of the field square from Ref. in the general case of the extra space we derived formulae 

for the VEV of the energy-momentum tensor in both regions on the right and on the left from 
the brane, given by expressions (|2.i;-ij) and (|2.2:-{j) respectively. Unlike to the case of purely AdS 
bulk, here the VEVs in addition to the distance from the brane depend also on the position of 
the brane in the bulk. In the limit when the AdS curvature radius tends to infinity we derive the 
formula for the vacuum energy-momentum tensor for parallel plates on the background spacetime 
with topology R^^^'^^ x S. In this limit for a homogeneous internal space |^-component of the 
brane induced part in the VEV of the energy-momentum tensor vanishes. Further we have 
investigated various limiting cases when the expression for the brane induced VEV is simplified. 
For the points on the brane the vacuum energy-momentum tensor diverges. The leading term 
in the corresponding asymptotic expansion near the brane is given by formula 1)2. 33() . Near 
the brane the total vacuum energy-momentum tensor is dominated by the brane induced part 
and has opposite signs for Dirichlet and non-Dirichlet boundary conditions. Near the brane 
^- and ^-components of this tensor have opposite signs in the regions y < a and y > a. For 
C, > Cdi (C < C-Di) the energy density is positive (negative) for Dirichlet boundary condition 
and is negative (positive) for non-Dirichlet boundary condition. For large distances from the 
brane in the region z > Za the contribution of a given mode along S with nonzero KK mass is 
exponentially suppressed by the factor e~'^^f^. For the zero mode the brane induced VEV near 
the AdS horizon behaves as z^^~'^'^. In the purely AdS bulk {D2 = 0) this VEV vanishes on 
the horizon for 1^ > 0. For an internal spaces with D2 > 2^ the VEV diverges on the horizon. 
The VEV integrated over the internal space (see Eq. (|2.26() ) vanishes on the AdS horizon for 
all values D2 due to the additional warp factor coming from the volume element. For the points 
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near the AdS boundary, the brane induced VEV vanishes as z '^'^'^ for diagonal components and 
as for the ^-component. For small values of the length scale for the internal space, the 

contribution of nonzero KK masses is exponentially suppressed and the main contribution into 
the brane induced energy-momentum tensor comes from the zero mode. In the opposite limit, 
when the length scale of the internal space is large, to the leading order the vacuum energy- 
momentum tensor reduces to the corresponding result for a brane in the bulk AdSo+i given in 
Ref. |17j . For strong gravitational fields corresponding to small values of the AdS curvature 
radius, the contribution from nonzero KK modes along S is exponentially suppressed by the 
factor e"^^''!^^^"!. For the zero KK mode the components of the brane induced vacuum energy- 
momentum tensor behave like k^^~^^e^^^°y exp[(L)i + 2v)kj:){y — a)] in the region y < a and like 
k^^~^^ e^^^^y eKp[2ukD{a — y)] in the region y > a. The corresponding quantities integrated over 
the internal space contain additional factor e~^^^^y coming from the volume element and are 
exponentially small in both regions. For fixed values of the other parameters, the brane induced 
VEV in the region z > Za vanishes as when the brane position tends to the AdS boundary. 
When the brane position tends to the AdS horizon, — > oo, for massive KK modes along S the 
VEV of the energy-momentum tensor in the region z < Za is suppressed by the factor e"^^"'^/'. 
For the zero mode in the same limit the suppression is power-law with respect to Za- 

The geometry of two branes we consider in sectional The VEVs in the region between the 
branes is presented in the form (|3.5() with separated boundary-free, single branes and interference 
parts. The latter is given by formula (|3.6|) and is finite everywhere including the points on the 
branes. The surface divergences are contained in the single brane parts only. We have explicitly 
checked that the both single brane and interference parts separately satisfy the continuity equa- 
tion and are traceless for a conformally coupled massless scalar. The possible trace anomalies 
are contained in the boundary- free parts. In the limit ^ we derive the corresponding 
results for two parallel Robin plates in the bulk R^^^'^^ x S. For small values of the length scale 
of the internal space corresponding to large KK masses, the interference part in the VEV of the 
energy-momentum tensor is estimated by formula (|3.9|) and is exponentially suppressed. For 
large interbrane distances compared with the AdS curvature radius kjy^ we have approximate 
formula (|.S.in|) for the nonzero KK modes along E. This limit is realized in braneworld scenarios 
for the solution of the hierarchy problem. The interference part vanishes as z'^'^ when the left 
brane tends to the AdS boundary. Under the condition z <^ Zb an additional suppression factor 
appears in the form {z/ z^)^^ for ^-component and in the form (z/zf,)^^'*'^"! for the other com- 
ponents, where ai = min(l, u). In a higher dimensional generalization of the Randall-Sundrum 
braneworld based on the bulk AdSoi+i x S with orbifolded y-direction the VEV of the energy- 
momentum tensor is obtained from the formulae in this paper with an additional factor 1/2. 
For untwisted scalar field the corresponding Robin coefficients are related to the surface mass 
parameters by formulae ()3.13|) . For twisted scalar Dirichlet boundary conditions are obtained. 
Note that in the present paper we consider the VEV of the bulk energy- momentum tensor. On 
manifolds with boundaries in addition to this part, the energy- momentum tensor contains a 
contribution located on on the boundary (for the expression of the surface energy-momentum 
tensor in the case of arbitrary bulk and boundary geometries see Ref. jH]). As it has been 
discussed in Refs. (22 |^ [2^1 ESl 123 EE], the surface part of the energy-momentum tensor is 
essential in considerations of the relation between local and global characteristics in the Casimir 
effect. The vacuum expectation value of the surface energy-momentum tensor for the geometry 
of two parallel branes in AdSo+i bulk is evaluated in Ref. [20] • In particular, it has been 
shown that for large distances between the branes the induced surface densities give rise to an 
exponentially suppressed cosmological constant on the brane. The investigation of the surface 
densities and induced cosmological constant on the branes for the bulk AdS^i+i x 5] will be 
reported in the forthcoming paper. 

The vacuum effective pressure acting on the branes is determined by the ^-component of 
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the energy-momentum tensor and can be separated into single brane and second brane induced 
parts. The first one is divergent and needs additional renormalization. The corresponding 
procedure lies on the same line as the evaluation of the surface energy-momentum tensor and 
will be discussed in the forthcoming paper. Here we concentrate on the second brane induced 
part which is finite for all nonzero interbrane distances and is not changed by the regularization 
and renormalization procedure. For the brane at z = zj this term is determined by formula 1)4. 
For Dirichlet scalar the corresponding vacuum forces are attractive for all interbrane distances. 
Taking the limit A:/) — > we obtained the result for the bulk x S, given by the right 

hand side of Eq. (|4.5j) . In this case, for a homogeneous internal space the interaction forces 
are the same even in the case of different Robin coefficients for separate branes. For the modes 
along E with large KK masses, the interaction forces are exponentially small. In particular, for 
sufficiently small length scale of the internal space this is the case for all nonzero KK modes and 
the main contribution to the interaction forces comes from the zero mode. For small interbrane 
distances, to the leading order the interaction forces are given by formula (|4.8|) . In this limit 
they are repulsive for Dirichlet boundary condition on one brane and non-Dirichlet boundary 
condition on the another and are attractive for other cases. For small interbrane distances the 
contribution of the interaction term dominates the single brane parts, and the same is the case 
for the total vacuum forces acting on the branes. When the right brane tends to the AdS horizon, 
Zb — > oo, the interaction force acting on the left brane vanishes as e~'^^f^^^ j z^^^"^ for the nonzero 
KK mode and like z'^^^~'^" for the zero mode. In the same limit the corresponding force acting 
on the right brane behaves as 2;^2+^i/2+ig-2A^2;j, £qj, ^.j^g nonzero KK mode and like z^'^~'^^ for 
the zero mode. In the limit when the left brane tends to the AdS boundary the contribution of 
a given KK mode into the vacuum interaction force vanishes as z^^'^^ and as z^ for the left and 
right branes, respectively. For small values of the AdS curvature radius corresponding to strong 
gravitational fields, for nonzero KK modes under the conditions XpZa ^ 1 and \(i{z\, — Za) 3> 1, 
the contribution to the interaction forces is suppressed by the factor e^^^^''^^''^^"). For the zero 
KK mode, the corresponding interaction forces integrated over the internal space behave as 
k^^~^^ exp[{Di6j + 2h')k£i{a — b)] for the brane at y = j and are exponentially small. In the 
model without the internal space we have similar behavior with v replaced by ui and for a scalar 
field with C < Cz)+Di+i the suppression is relatively weaker. 

As an example of application of the general results, in section |S1 we consider a special case 
of the internal space with T, = S^. By using the Abel-Plana summation formula, the boundary- 
free part of the vacuum energy-momentum tensor is presented in the form 1)5. 2() , where the first 
term on the right is the corresponding tensor in AdS^+i bulk without boundaries. The second 
term is induced by the compactness of the X direction and behaves as z^^'^'^ for z <^ L and 
like z^~^^ for z ^ L. For minimally and conformally coupled scalar the corresponding energy 
density is positive in the first case and negative in the second case and has a maximum at some 
intermediate value for z/L. The total energy- momentum tensor is dominated by the first term on 
the right of Eq. (|5.2() near the AdS boundary and by the second term near the horizon. For the 
points near the branes the brane induced parts dominate. In the case of internal space S = 5^ 
the latter are obtained from the general formulae by replacements 1)5. 9() . We also comment on 
the generalization to the case of the internal space 
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